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ABSTRACT 

The turnover of interstellar gas on ~ 10 9 yr timescales argues for the continu- 
ous operation of a galactic dynamo. The conductivity of interstellar gas is so high 
that the dynamo must be "fast" - i.e. the magnetic field must be amplified at a 
rate nearly independent of the magnetic diffusivity. Yet, all the fast dynamos so 
far known - and all direct numerical simulations of interstellar dynamos - yield 
magnetic power spectra that peak at the resistive scale, while galactic magnetic 
fields have substantial power on large scales. In this paper we show that in 
weakly ionized gas the limiting scale may be the ion-neutral decoupling scale, 
which although still small is many orders of magnitude larger than the resistive 
scale. 

Subject headings: dynamos — MHD — methods: numerical — ISM: magnetic fields 

1. Introduction 

Despite many theoretical and observational advances, our understanding of galactic 
magnetic fields is still incomplete. Although there is evidence for magnetic fields in young 
galaxies (Kulsrud & Zweibel 2008), it is likely that dynamo processes still operate continu- 
ously in galaxies today. Perhaps the most compelling argument for ongoing dynamo activity 
is that the turnover time of interstellar gas due to loss and replacement of material - ~ 10 9 
yr in the case of the Milky Way - is much less than the ages of galaxies. As interstellar gas 
is added - whether by infall from the IGM or through stellar mass loss - its magnetic field 
must be brought to the galactic value so as to maintain the field in a steady state. 
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In contrast to the primordial situation at early times, when a magnetic field had to be 
built up from nearly zero, it is likely that all the material on which galactic dynamos now 
act is at least somewhat magnetized. Under these conditions, the primary functions of the 
dynamo are to increase magnetic energy to the observed value, to generate and maintain 
a component of the field which is coherent over at least several kpc, as observed, and to 
transport magnetic flux into the gas which has newly arrived, which especially if it is in- 
tergalactic in origin or was shed by stars, may be under-magnetized. These processes are 
also necessary in theories such as that of Rees (1987), in which galactic magnetic fields are 
seeded by random fields injected by many small scale sources. 

Dynamos are considered "fast" if the rate at which the field is amplified tends to a 
finite limit as the magnetic diffusivity r\ approaches zero (Childress & Gilbert 1995). It may 
seem paradoxical that magnetic diffusivity, which removes energy from the magnetic field, is 
required for dynamos at all. But without diffusion the magnetic topology is fixed, and this 
places strong constraints on magnetic field amplification (Moffat 1978). The Ohmic diffusion 
time in galaxies is so much longer than any dynamical time that it seems galactic dynamos 
must be fast. 

Fast dynamo theory has followed two approaches. One involves computation and anal- 
ysis of the action of prescribed flows on magnetic fields under nearly non-diffusive conditions 
(e.g. Galloway & Proctor 1992, Ott 1998, Gilbert 2002, Couvoisier, Hughes, & Tobias 2006). 
These so-called kinematic studies have demonstrated the role of chaotic flow in magnetic field 
amplification and elucidated important relationships between the topological properties of 
the flows and their dynamo properties (Klapper & Young 1995). The modification of some 
of these flows by magnetic forces has also been considered, and has provided some insight 
into how dynamos might saturate (Cattaneo, Hughes, & Kim 1996, Tanner & Hughes 2003, 
Cameron & Galloway 2005). 

The second approach to fast dynamos involves analysis and direct numerical simula- 
tions of driven magnetohydrodynamic (MHD) turbulence, again with the lowest magnetic 
diffusivity A possible. In particular, A is chosen to be less than the viscous diffusivity v (i.e. 
the large Prandtl number case). These models are self-consistent in the sense that magnetic 
forces are fully included in the dynamics, and theory and simulation are in good agreement 
on how the field evolves (Schekochihin et al. 2004). 

Both approaches yield magnetic fields which are dominated by structure on the resistive 
scale, which is far smaller than any other characteristic scale in the interstellar medium. This 
follows from the fact that in the absence of magnetic diffusion, any divergence free flow which 
amplifies the field lengthens the field lines in the same proportion. Therefore, as the field 
is amplified it becomes folded or tangled. The large random field is in stark contrast to the 
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observed structure of galactic magnetic fields, which display considerable long range order, 
and thus presents a serious challenge to dynamo theory. In particular, the magnetic field is 
amplified on scales far below the minimum velocity scale, whether this scale is prescribed 
(the kinematic case) or created by viscous effects (the full turbulent simulation). 

Most of the mass, and a considerable part of the volume of the interstellar medium, 
is weakly ionized, with the plasma and neutral fluids coupled by collisions and by ioniza- 
tion and recombination. On long lengthscales and timescales the medium behaves like a 
single conducting fluid, but on short lengthscales and timescales, the two fluids decouple. 
For present interstellar medium parameters, the magnetic field is in approximate energy 
equipartition with the bulk fluid and dominates the plasma component. The question there- 
fore arises whether a dynamo in a weakly ionized medium could be quenched on scales below 
the plasma-neutral decoupling scale while continuing to amplify the field on larger scales. 
Based on analysis and computations presented in this paper, the answer seems to be that it 
can. 

This is far from the first study of the effects of partial ionization on dynamos. There 
have been a number of studies of large scale dynamos in weakly ionized gases which show 
that the nonlinearity introduced by ion-neutral drifts has important consequences for mean 
field theory (Zweibel 1988, Proctor & Zweibel 1992) and for the general problem of how large 
scale fields are amplified in a medium with small scale turbulence (Kulsrud & Anderson 1992, 
Subramanian 1997, 1999, Brandenburg & Subramanian 2000). These studies all used the 
strong coupling approximation (Shu 1983), according to which the plasma drift relative to 
the neutrals is found by balancing the Lorentz force by ion-neutral drag. This approximation 
breaks down at small scales. In this work, we treat the plasma and neutrals as two frictionally 
coupled fluids, which is more appropriate at the small scales, and also more general. 

In order to see how ion-neutral friction might affect the dynamo, consider the equation 
for magnetic energy, which in a domain V with periodic or infinitely distant boundaries can 
be written as 



where Vj is the plasma velocity and A is the Ohmic diffusivity. The first term in the integral 
on the right hand side represents the work down by the field on the plasma; the second term 
is resistive dissipation. 

The plasma velocity can always be expressed in terms of the neutral velocity v n and 
the drift velocity Vd = v* — v n ; V; = v„ + v D . If the ionization fraction is low, v„, in turn, 
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is very nearly the center of mass velocity v. In the strong coupling approximation, vp is 

_ (V x B) x B 
v £> - — : > K z ) 

where 7 = (&v) in / (rrii + m n ) is the ion-neutral friction coefficient. Replacing Vj by v + v D 
in equation (JTJ and using equation (J2J), we see that when the strong coupling approximation 
holds, the magnetic energy equation is 
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Equation ([3]) shows, as expected, that ion- neutral friction is a dissipative effect. The work 
term, however, can have either sign, and it is possible that in the presence of ambipolar 
drift, the magnetic field and flow are modified so as to increase the rate at which energy 
flows to the field. Equation (pQ), on the other hand, is valid whether the strong coupling 
approximation holds or not, is closer to the magnetic energy equation for a plasma, and does 
not attempt to split the inductive and dissipative effects. 

In §[2] we introduce the basic setup. In §|3]we show how plasma- neutral friction can allow 
the field to saturate at the small scales while still being amplified at the large scales. Section 
H]is a summary. 



2. Formulation 

2.1. Important timescales and lengthscales 

In a turbulent cascade, the velocity at scale I, vi, typically decreases with decreasing I, 
but slowly enough that the turnover time 77 = l/vi, decreases with I as well. For example, 
in Kolmogorov turbulence, or in the magnetohydrodynamic turbulence model of Sridhar 
& Goldreich (1994) and Goldreich & Sridhar (1995, 1997), v\ is related to / and the scale 
Id at which the turbulence is driven by vi ~ Vd{l/ld) 1 ^ ■ The turnover time 77 is then 
(ld/vd)(l/ld) 2 ^ 3 - The cascade terminates at the scale Ik at which the dissipation rate becomes 
faster than the turnover rate. 

In diffuse, weakly ionized interstellar gas, tk, the turnover time at Ik, is comparable 
to the neutral-ion collision time r n ,, and much longer than the ion-neutral collision time 
Ti n . Only for motions on timescales less than Ti n is friction with the neutrals unimportant 
for plasma dynamics. Resistive effects are typically unimportant for motions with turnover 
times r in , but become important on much smaller scales. The ion viscous scale is probably 
terminates the ion flow on scales much larger than the resistive scale (the suppression of 
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viscous stresses perpendicular to B even for a weak magnetic field makes the actual value 
difficult to assess), and the 

The computational resources available to us preclude modeling all these scales in a 
realistic manner. Our main focus is on separating the ion-neutral decoupling scale and the 
resistive scale. As a consequence, instead of assigning the neutral motions a third, still 
larger scale, we allow the neutral flow to take place on the ion-neutral decoupling scale. 
Furthermore, we rely entirely on numerical viscosity, which allows the flow in the ions to 
extend below the resistive scale. In reality, ion viscosity probably terminates the ion flow on 
scales much larger than the resistive scale (the suppression of viscous stresses perpendicular 
to B even for a weak magnetic field makes the actual value difficult to assess), but again we 
are unable to achieve true scale separation (or implement anisotropic viscosity). 

The consequences of these aspects of our calculations are discussed in the next two 
sections. 



We assume that the neutral velocity v n is of a form which would give fast dynamo 
action if it were the plasma velocity Vj. On timescales much longer than the ion-neutral 
collision time (pnlin) -1 , the plasma should move with the neutrals - Vj ~ v ra - and the 
magnetic field should grow at the fast dynamo rate, at least while Lorentz forces are unim- 
portant. On short lengthscales, inertial and Lorentz forces compete with ion- neutral friction 
to drive Vj away from v n . We investigate the dynamo properties of the resulting plasma 
flow. The setup is somewhat similar to the two-fluid study of Kim (1997), although that 
work focused on diffusion of a large scale field, and the computations were two dimensional 
and incompressible. 

We simplify the problem by assuming the ionization and recombination timescales are 
much shorter than the advective timescales, so that ionization equilibrium holds. This 
assumption allows us to bypass the plasma continuity equation, but it is unrealistic in the 
diffuse gas which is the primary application for this study [although it is much better for 
molecular gas (see Table 3 of Heitsch & Zweibel 2003a)]. The neutral velocity we choose is 
divergence free, so by taking the neutral density p n to be initially uniform, we force it, and 
the plasma density p p ~ pi, to remain so. Therefore, there are no thermal pressure gradient 
forces in the problem. Under these conditions, the momentum equation for the plasma is 



2.2. Case study 
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where v n is to be prescribed. Equation (j3J) is to be solved together with the magnetic 
induction equation 

<9B 

_ = V x ( v< x B) + AV 2 B. (5) 

For v n we choose the 2.5D flow shown by Galloway & Proctor (1992) to be a fast 
dynamo. This flow, which we will refer to as v G p, can be written as 

v n = x2Tck GP 4) + y— - z—, (6) 

where the stream function ip is 

--—r- — ( sin [2%(k G pz + esin27ra;t)] 
2 2nkcp V 

+ cos [2n(kGpy + e cos 2itut)\ J , (7) 

and e is a constant. This flow has periodic, cellular structure on a single spatial scale, k^ 1 , 
and temporal structure at frequency to and all of its harmonics. This can readily be seen by 
rewriting equation ([7j) in the form 



oo 

ip = ip n ( sin [2Tc(k GP z + nut)] 



where 



n=0 

+ cos[2n(k GP y + n( y ut+ -)]), (8) 

and the J n are ordinary Bessell functions of the first kind. However, bearing in mind that 
Jn{%) ~ (x/2) n /n\ for x < n, and that x = n/2 in our problem, we see that only the first 
few harmonics of the series are important in the flow. The cell pattern drifts back and forth 
with temporal frequency to and amplitude e. The flow has chaotic regions along the cell 
edges. In this paper, we choose units of length and time such that k G p — Vq — u = 1. The 
eddy turnover time and the period over which the pattern oscillates are then of the same 
order. We set e = 0.25, near the value at which chaos is maximized (Brummell, Cattaneo, 
& Tobias (2001)). The ranges of y and z are —1 < y < 1, — 1 < z < 1, so only one GP cell 
fits into the domain. We choose A such that Xk GP , the resistive decay rate at the GP scale, 
is between 10~ 2 and 10~ 4 , and the friction coefficient, 7i„p n , is between 0.03 and 0.60. We 
set pi — 1. The magnetic field is initialized by choosing a vector potential which is a Fourier 
series in modes of the computational domain with random amplitudes and random phases 
selected from uniform distributions. 
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The numerical scheme, Proteus, is based on the conservative gas-kinetic flux splitting 
method introduced by Xu (1999) and Tang & Xu (2000), and is second order in time and 
space. Resistivity is included through dissipative fluxes (Heitsch et al. 2004, 2007), but 
viscosity is not (there is a small amount of numerical viscosity). The frictional force is 
implemented through a drag term (Heitsch et al. 2004). 



3. Results 



The complete set of models we considered is summarized in Table 1. We first verified 
that our computational scheme recovers the previously known dynamo behavior of the GP 
flow, and assessed numerical convergence. The left panel of Figure [1] shows the magnetic 
energy against time for three resolutions and five diffusivities. Once the short wavelength 
magnetic power present in the random initial conditions has decayed, the magnetic energy 
grows at an exponential rate. The growth rates are plotted in the right panel of Figure [TJ 
With our choice of units, R m = A -1 . While some flattening of the curves is apparent, the 
growth rate has evidently not converged to the R m — > oo limit. Nevertheless, the evident 
lack of convergence at the largest _R m and highest numerical resolution practical for us shows 
that it would not be meaningful to set A to any value smaller than what we used here. This 
series of models will be referred to as KI (for model keys, see Table [T]). All the models 
discussed subsequently have R m = 10 4 . 
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Fig. 1. — Left: Natural logarithm of the normalized magnetic energy against time for all 
kinematic models (series KI), at linear resolutions of N = 64 (thin lines), 128 (medium lines), 
and 256 (thick lines). The higher A, the stronger the initial decay of the (noise-like) seed 
field. Right: (Exponential) growth rates derived from fits to the curves in the left panel, for 
all resolutions and magnetic Reynolds numbers. 
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Table 1. Key to model names 



Name 


Mnemonic 


Physics 


KI 


Kinematic, Ions 


single-fluid, GP flow in ions 


KN 


Kinematic, Neutrals 


two-fluid, GP flow in neutrals 


DL 


Dynamic, Lorentz 


two-fluid, GP flow in neutrals, Lorentz force 


DR 


Dynamic, Reynolds 


two-fluid, GP flow in neutrals, Reynolds Stress 


DA 


Dynamic, All 


DL combined with DR 
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We next set v n = v^p and solved equation with the Lorentz force and Reynolds 
stress terms dropped. These models, which we denote by KN, demonstrate the interplay 
between inertia and friction in determining the ion flow. Using equation (jSj), equation (j3j) 
can then be reduced to a series of equations for the jth Fourier harmonic of 



which can be solved analytically. For the parameters chosen here, p n ji n <C joo for j > 0, and 
is nearly 7r/2 out of phase with the frictional forcing term, and reduced in amplitude by a 
factor of Pnjin/ju. Thus, inertial effects force the plasma away from the GP flow (this is an 
artifact of the overlap in our models between the turnover rate of the neutral eddies and the 
ion-neutral collision rate; in the interstellar medium the latter is much larger). However, even 
though the plasma flow departs significantly from the GP flow, it does amplify the magnetic 
field. This is shown in Figure EJ The amplification rate decreases with decreasing p„7j n . 
This is caused primarily by the reduced amplitude of the ion flow, but also by its structure, 
which as Pnlm decreases departs more and more from w GP . The dissipative nature of ion- 
neutral friction may also be playing a role here; see equation - although equation (|2J) is 
not completely accurate in this case. 

It can be seen from equation (jSJ) that magnetic fluctuations at wave number k interacting 
with flow at the GP wave number k G p are coupled to fields at wave number k ± kcp- This 
sets off a magnetic cascade in wavenumber, and therefore, both the GP flow and the modified 
flow present in the KI models drive magnetic field at all spatial scales. 




(10) 
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Fig. 2. — Natural logarithm of the normalized magnetic energy against time for models at 
linear resolution N = 128. See Table [T] for keys to the model names. Back-reaction via 
Lorentz- forces let the field amplification saturate (models KN vs DL), while the inclusion of 
the Reynolds stress suppresses the dynamo (models KN vs DR and DA). Note that there 
is no difference in the energy evolution between models DR and DA because the magnetic 
field simply stays too weak for the Lorentz force to be of importance. The dashed horizontal 
line denotes energy equipartition between magnetic energy and the nominal kinetic energy 
in the GP flow. 
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The Reynolds stress Vj ■ Vvj, because of its nonlinearity, drives higher spatial harmonics 
of the GP wave number. In our parameter regime, it is of the same order as the inertial 
term, at least if t>j is of the same magnitude as vqp- Because in our calculation the re- 
sistivity exceeds the (numerical) viscosity, the ion flow extends to smaller scales than the 
magnetic field. It is known that dynamo activity is often suppressed in these so-called low 
magnetic Prandtl number situations unless the magnetic Reynolds number is large (Boldyrev 
& Cattaneo 2004, Schekochihin et al 2007). Because the ion density remains constant, the 
compressibility of the medium is exaggerated; suppression of the dynamo has also been as- 
sociated with compressibility effects (Haugen, Brandenburg, & Mee, 2004). For both these 
reasons, it is perhaps unsurprising that including the Reynolds stress obliterates the dynamo; 
see the DR models in Figure [2j 

We also considered the effect of dropping the Reynolds stress term but including the 
Lorentz force (models DL). Because it is nonlinear in B, and because magnetic power is driven 
at all scales, the Lorentz force too drives higher spatial harmonics in the flow. The feedback 
from the magnetic field eventually saturates the dynamo (Fig. [2]). At the time of saturation, 
the volume integrated plasma kinetic energy is about 100 times larger than the magnetic 
energy. However, magnetic forces along the cell walls, where the field is concentrated (see 
Figure [3]), are locally strong, and inhibit the exponentially fast stretching needed for further 
amplification. This is consistent with the results of Cattaneo, Hughes, & Kim (1996), who 
showed that magnetic forces suppress chaos - characterized by exponential stretching - in 
the Galloway-Proctor flow even when the field is too weak to modify the flow kinetic energy. 

Models with the most complete physics, namely both Reynolds stress and Lorentz forces 
in addition to inertia and ion- neutral friction, fail to be dynamos (Fig. [2j models DA). This 
is to be expected, since the initial magnetic field is too weak to affect the dynamics, and 
nonlinear advection alone suppresses the dynamo. 

It is revealing to look at the structure of the magnetic field as well as its growth rate. 
Figure [3] shows the logarithm of the magnetic energy density from three views (along the x, 
y, and z axes) for the KI, KN, DL, DR, and DA models, at t — 20. The top set of panels 
represent the canonical kinematic GP dynamo (KI). The dark lines, where the magnetic field 
is very small, demonstrate that it is folded on scales much smaller than the flow scale, as 
is characteristic of dynamos with large P m . The second row of panels shows the field when 
the plasma flow is determined by friction and inertia, but not Lorentz forces or Reynolds 
stresses (KN). It is weaker than in models KI, but still spatially intermittent and folded. If 
these folds were true magnetic nulls, they could be sites of fast magnetic reconnection, but 
even a small amount of magnetic shear is enough to quench this process (Heitsch & Zweibel 
2003a,b). 




Fig. 3. — Logarithm of magnetic energy density at midplane, in planes (y, z) (left), (x,z) 
(center) and (x, y) (right). Models from top to bottom: KI, KN, DL, DR and DA. For model 
names, see Table CD 



-14- 



More dramatic changes come with the inclusion of the Lorentz force, but not the 
Reynolds stress (third row of Figure 3, models DL). Although Figure [2] shows that the 
magnetic field is still some way from saturation, and the magnetic energy is indistinguish- 
able from that of the kinematic models, the spatial intermittency and amount of folding are 
markedly reduced. 

The final two rows, models DR and DA, show how the small scale flow generated by 
nonlinear advection affects the magnetic field. The ordered orientation associated with the 
strong shear in the GP flow is essentially gone. Since there is little dynamo action, magnetic 
feedback is unimportant, and the magnetic fields in the two cases are virtually identical. 
Although there has been little amplification of the field as a whole, there are thin filaments 
where it is locally strong. These structures are associated with strong compression (negative 
divergence in the ion velocity) and with strong vorticity. 

Magnetic energy spectra for the KI, KN, and DL models at a linear resolution of iV = 256 
are shown in Figure HI Although there is a slight reduction of power at the smallest scales 
it is not as important as the overall diminished amplitude and it is clear that power spectra 
without phase information give an incomplete picture of the structure of the magnetic field. 
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Fig. 4. — Power spectra of the magnetic energy for models KI, KN and DL at a linear 
resolution of iV = 256, and at system time t = 16, i.e. at a point when the energy evolution 
of models KN and DL has separated due to saturation (see Fig. [2]). 
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4. Discussion and Summary 

Galactic magnetic fields show considerable long range order, while computations of small 
scale turbulent dynamos in highly conducting plasmas have shown instead that magnetic 
energy is concentrated at the resistive scale. This suggests that there must be some feature 
of the interstellar dynamo which suppresses growth of the small scale field, and that there 
exist mechanisms for transferring energy to scales larger than the turbulent injection scale. 

In this paper we investigated the role of ion-neutral friction in eliminating the smallest 
scales. We solved a simple model problem in which the neutral flow was known to be a 
dynamo, and was coupled to the plasma through frictional drag. The plasma flow was 
determined by friction, inertia, and Lorentz forces. Ionization equilibrium was assumed. 
We used a spatially constant Ohmic resistivity, and we relied on viscous dissipation at the 
resolution scale. 

This study is, as far as we know, the first dynamo model in which the plasma flow is 
calculated from the full equation of motion, rather than assuming the strong coupling ap- 
proximation (Shu 1983). The strong coupling approximation breaks down at small scales, 
and leads to an equation for magnetic energy growth which differs somewhat from the stan- 
dard equation in a plasma [compare equations (pQ) and ([31]. However, our models differ from 
the interstellar medium in several significant respects, due to the prohibitively large com- 
putational resources required for a realistic simulation. In the models, the neutral forcing 
occurs on timescales comparable to the ion-neutral collision time, whereas in the diffuse, 
weakly ionized interstellar medium the shortest timescale in neutral turbulence is probably 
2-3 orders of magnitude longer. This introduces an artificially large difference between Vj 
and v„. By taking the magnetic diffusivity larger than the viscous diffusivity (small magnetic 
Prandtl number P m ), we reverse the interstellar ordering (P m ^> 1). By assuming ionization 
equilibrium, an isothermal equation of state, and a uniform density of neutrals, we preclude 
the possibility of plasma pressure gradients. And finally, the resistive scale is only about two 
orders of magnitude less than the dynamical scale, instead of 5-6, as in the ISM. 

Our work does not address the growth of the field at the largest scales. Whether 
ambipolar diffusion acting on small scale turbulence can drive a large scale field, either by 
quasilinear (Zweibel 1988, Proctor Sz Zweibel 1992) or nonlinear (Subrahmanian 1997, 1999, 
Brandenburg & Subrahmanian 2000) effects cannot be probed by our study due to the large 
range of scales that would be required, and also, perhaps, because the neutral velocity is 
already a dynamo. 

Our main results are as follows. First, friction successfully imprints a dynamo flow on 
the ions, despite the importance of inertial terms [models KN; Fig. (j2J)]. That is, computing 
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Vj from equation ffTUj) leads to a flow which amplifies a weak magnetic field. Second, adding 
nonlinear advection (vj ■ Vv») to equation ffTUj) results in an ion flow which does not amplify 
the magnetic field, except in thin regions with large vorticity or negative divergence [models 
DR; Fig. Q]. This may be a manifestation of the difficulty of driving a dynamo when 
P m <C 1, or in a highly compressible medium. It would not necessarily carry over to the 
diffuse ISM, in which P m is larg^] and the recombination time is longer than the eddy 
turnover time. Third, when the Lorentz force is added, but the Reynolds stress is dropped, 
the dynamo saturates when the magnetic energy is about two orders of magnitude less than 
the plasma kinetic energy [models DL; Fig. (j2J)]. According to conventional wisdom, the 
magnetic energy saturates only close to equipartition; the difference here is due to the spatial 
intermittency of the magnetic field, which allows Lorentz forces to become strong just where 
the field is being amplified. Finally, adding the Lorentz force suppresses small scale structure 
(Fig. W)- 

The saturation reported here for models without nonlinear advection (models DL) has 
some features in common with the study by Tanner & Hughes (2003), although these authors 
considered fully ionized systems. They studied a superposition of the GP flow (which they 
call the CP, or Circularly Polarized flow) and a steady flow known to be a slow dynamo, 
restricted the dynamics by omitting the Reynolds stress term and averaged the Lorentz force 
over what in our case would be the x direction. They found that when the GP flow dominates, 
saturation occurs through suppression of exponential field line stretching by Lorentz forces, 
rather than by enhanced dissipation brought about by an increase in small scale structure. 

Taken at face value, our computations suggest that in weakly ionized gases, the efficiency 
of dynamos is reduced at scales below the neutral-ion coupling length. The main point is 
that saturation of the field at small scales is determined by the kinetic energy in the plasma 
at these scales rather than in the bulk medium. The inertia of the neutrals contributes, and 
presumably increases the saturation level, only for motions with turnover time t\ greater than 
the neutral-ion collision time r n j = (pi7i„) _1 . If we take rii ~ 2 x 10~ 3 cm -3 , representative 
of both cold and warm neutral interstellar gas, then t„; ~ 3 x 10 11 s. Taking the expressions 
for Kolmogorov turbulence given in §2.1, and assuming the driving scale and velocity are 10 
pc and 10 km s -1 , respectively, we find Ti ~ r ni at / ~ 10~ 2 pc. The equipartition magnetic 
field strength at this scale relative to the ions is only ~ .02fiG, or about 1% of the Galactic 
value. Thus, saturation of the dynamo below the neutral-ion decoupling scale could prevent 
amplification of the field at very small scales. 



lr The viscosity tensor in interstellar plasma is highly anisotropic due to the large ratio of the ion gy- 
rofrequency to the Coulomb collision rate. Therefore, viscous transport across the magnetic field is strongly 
suppressed (Braginskii 1965), making the actual value of P m somewhat unclear. 
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At early times, when the magnetic field might have been several orders of magnitude 
or more weaker than it is today, the decoupling effect would have been less important. Fur- 
thermore, the interstellar medium was probably almost fully ionized until metallic coolants 
had mixed into it. Thus, the mechanisms addressed here are probably more relevant to the 
maintenance of galactic magnetic fields than to their ultimate origin. 

We are happy to acknowledge supports by NSF Grants AST-0507367 and PHY-0215581 
to the University of Wisconsin. Parts of the computations were performed at the National 
Center for Supercomputing Applications under projects DAC AST-040026 and MRAC AST- 
060031. We are grateful to an anonymous referee for useful suggestions. 
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